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In QCD with two light quarks with apphcation of Bogolubov quasi-averages ap- 
proach a possibihty of spontaneous generation of an effective interaction, leading 
to the Nambu - Jona-Lazinio model, is studied. Compensation equations for 
form-factor of the interaction is shown to have the non-trivial solution leading to 
theory with two parameters: average low-energy value of as and dimensional pa- 
rameter /tt- All other parameters: the current and the constituent quark masses, 
the quark condensate, mass of vr meson, mass of a meson and its width are 
expressed in terms of the two initial parameters in satisfactory correspondence 
to experimental phenomenology. The results being obtained allow to state an 
applicability of the approach in the low-energy hadron physics and promising 
possibilities of its applications to other problems. 



1 Introduction 

In work [I] based on quasi-averages approach by N.N. Bogolubov |21ini a method is proposed 
aimed to obtain an effective interaction in a renormalizable quantum field theory. In par- 
ticular such an effective interaction is necessary for construction of the well-known Nambu 
- Jona-Lazinio model |3], which describes the low-energy physics of strong interactions. In 
view of this it is of interest to apply the developed method [T] to a study of possibility 
of generation of effective four-fermion interaction, which is intrinsic to the Nambu - Jona- 
Lazinio model. Following the results of work [T] one is to expect the essential contraction of 
number of initial parameters of the model. 

An application of the method [I] to a rather extensively studied low-energy region of 
hadron physics may show to what extent the use of the method is justified. In the present 
work the first non-perturbative approximation of the method will be developed in application 
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to the problem of a spontaneous generation of Nambu - Jona-Lazinio Lagrangian in the 
genuine strong interaction theory, namely in QCD with doublet of light quarks. 

2 Compensation equation for effective form-factor 

Now we start with QCD Lagrangian with two light quarks {u and d) with number of colours 
= 3 

2 

k=l 

where we use the standard QCD notations. 

In accordance to the Bogolubov approach, application of which to such problems being 
described in details in work Ij, we look for a non-trivial solution of a compensation equation, 
which is formulated on the basis of the Bogolubov procedure "add - subtract". Namely let 
us write down the initial expression in the following form 

L = Lq + Lint ] 

- - \ G / - - - - \ 

Unt = gsi^i,t^A;ij - i - F^Vi^oV) - y ' (v^As^V^As^ - ^^^v^) - 

Q / _ _ _ \ 

- -J- ■ yP'T''l^ip ipT^l^i' + ^r^757^V'?/'r%57^^ j . (3) 

Here ip is the isotopic doublet, colour summation is performed inside of each fermion bilinear 
combination, Fq^^, = d^A„ — diyA^, and notation Gi ■ ipiljipilj corresponds to non-local vertex 
in the momentum space 

t (27r)^ Gi Fi {pl,p2, p3, pA) 6{pl + p2 + p3 + p4) ; (4) 

where Fi{pl,p2,p3,pA) is a form-factor and pi, p2, p3, p4 are incoming momenta. In 
the same way we define vertices, containing Dirac and isotopic matrices. We comment the 
composition of the vector sector, being proportional to G2, in what follows. 

Let us consider expression (j21) as the new free Lagrangian Lq, whereas expression (jH)) 
as the new interaction Lagrangian Lint- Then compensation conditions (see again jT]) will 
consist in demand of full connected four-fermion vertices, following from Lagrangian Lq, to 
be zero. This demand gives a set of non-linear equations for form-factors F^. 

These equation according to terminology of works j2l 0] are called compensation equa- 
tions. In a study of these equations it is always evident the existence of a perturbative trivial 
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solution (in our case Gj = 0), but, in general, a non-perturbative non-trivial solution may 
also exist. Just the quest of a non-trivial solution inspires the main interest in such problems. 
One can not succeed in finding an exact non-trivial solution in a realistic theory, therefore 
the goal of a study is a quest of an adequate approach, the first non-perturbative approxi- 
mation of which describes the main features of the problem. Improvement of a precision of 
results is to be achieved by corrections to the initial first approximation. 

Thus our task is to formulate the first approximation. Here the experience acquired in 
the course of performing of work [Tj could be helpful. Now in view of obtaining the first 
approximation we assume: 

1) In compensation equation we restrict ourselves by terms with loop numbers 0, 1, 2. 

2) In expressions thus obtained we perform a procedure of linearizing, which leads to linear 
integral equations. It means that in loop terms only one vertex contains the form-factor, 
being defined above, while other vertices are considered to be point-like. In diagram form 
equation for form-factor Fi is presented in Fig 1. 

3) While evaluating diagrams with point-like vertices diverging integrals appear. Bearing in 
mind that as a result of the study we obtain form-factors decreasing at momenta infinity, we 
introduce a cut-off A in the diverging integrals. It will be shown that results do not depend 
on the value of this cut-off. 

4) We can obtain analytic expressions for massless quarks only. We take masses into account 
by introducing the cut-off in the lower limit of integration by momentum squared at a 
value, which equals to a value of the corresponding propagator denominator at g = 0. Namely 
with the denominator (g^ -|- m?) the cut-off parameter equals m?. In doing this we keep at 
nominators only the leading terms in m expansions because taking into account of the next 
terms evidently means supererogation of accuracy. 

5) We shall take into account at most the first two terms of the expansion. 

Let us formulate compensation equations in this approximation. For free Lagrangian 
Lq full connected four-fermion vertices are to vanish. One can succeed in obtaining ana- 
lytic solutions for the following set of momentum variables (see Fig. 1): left-hand legs have 
momenta p and — p, and right-hand legs have zero momenta. In particular this kinematics 
suits for description of zero-mass bound states. Under some assumptions solutions obtained 
under these conditions may be generalized to momentum set (p, —p,q, — g). The construc- 
tion of expressions with an arbitrary set of momenta is the problem for the subsequent 
appr oximat ions . 

Now following the rules being stated above we obtain the following equation for form- 
factor Fi{p) in scalar colour singlet channel 



GiF,{p') 




27r2 V AN 



GiN 



2n^ 




+ 



3Gi G2 
87r2 



2p^) 32 7r6 



/i2 X _ {Gl + QGiG2)N 
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X 



j(2A +ip-q) log^^^ -ip-q) -^(^j ; (5) 



Here integration is performed in the four-dimensional Euclidean momentum space, fi = m^. 
One-loop terms contains terms proportional to N and 1 while two-loop terms correspond 
to A'"^ and N . The leading terms are the same for scalar and pseudo-scalar cases. We 
begins the study with the scalar channel, because it defines chiral symmetry breaking effect. 
Equation (0) evidently has trivial solution Gi = 0. Bearing in mind our goal to look for 
non-trivial solutions we divide the equation by Gi. As a result of four-dimensional angle 
integration (see Appendix) we have 

3G2^^,2,^i_^ 3 ^2 {Gl + QGiG2)N 

^ X 1 ;; X 



Fi{x) = A + —^(2k^ + x\og- 



A2 2 2x) 32 TT^ 

3 



X ( ^ / (y^- 3/i')i^i(2/) dy + - yFiiy) dy + \ogx yFi{y) dy + 



+x\ogxj Fi{y)dy + j y log y Fi{y) dy + x j (log y + Fi{y) dy + 

POO POO 

- log A^U yFiiv) dy + X F^{y) dy 

GlNA^ 1 GiAT 1 X , ^ , , ^ > 



(6) 



X = v^] y = 



In view of taking into account of quark mass the lower limit of the momentum integration 
/i = mQ is introduced. Equation ® by a sequential six-fold differentiation reduces to the 
following differential equation 



iGl + 6GiG2)N 



/3 



16 vr^ 



with boundary conditions to be formulated below. 

Equation ((Tj) reduces to Meijer equation 0, IZ|- Namely with the simple substitution we 
have 

4 - ") (4 - <•) ^ i (4 4) (4 - ^) (4 - - ^) w = " ^ w 

px^ 1 - VI - 64m , 1 + VI - 64m pfi^ 

= — — : a = : o = : u = . 

26 ' 4 ' 4 ' 64 
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Boundary conditions for equation (jHJ are formulated in the same way as in works PP , j3] . At 
first we have to choose solutions decreasing at infinity, that is combination of the following 
three solutions 

F,{x) = C^Gtl[z\l, i i 0, a, 6) + C^Gt^^zll, ^, 6, a, ^, 0, ) + 

+ C3G^°(.|l,0,6,a, 1 i); ^ = ^- (9) 



where 



is the Meijer function [H], j7j with sets of upper indices Oj and of lower ones hj. In case only 
one line of parameters is written this means the presence of lower indices only, n and p in 
the case being equal to zero. 

Constants Cj are defined by conditions 

r2 



7r2 2 



Fiiy)dy = 0; 

POO 

/ yF,iy)dy = 0; (10) 

/>oo 

These conditions and condition A = as well provide cancellation of all terms in Eq. (jH)) 
being proportional to and log A^. Thus the result does not depend on a value of parameter 
A. By solving linear set (fTUj) . in which solution ^ is substituted, we obtain the unique 
solution. Value of parameter uq, which is connected with current quark mass, and ratio of 
two constants Gi we obtain from conditions = 1 and 

^ GiNA^ / 1 GiN 1 \ r \ 

= fl + ^l^fl 6G^(4iV + 2L^ 

V^ANJ 27r2 I (Gi + 6G2)(4iV+l)y ' 

this gives for = 3 with the account of the first of conditions (fTUI) 

no = 1.6-10-«; d = ^ . (12) 

The form-factor now reads as Q with 

Ci = 0.28322 ; C2 = - 3.655 ■ 10~^ ; C3 = - 7.794 ■ 10^^ ; (13) 
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-^i('^^o) = 1 Fi{z) decreases with z increasing. It is important, that the solution exists 
only for positive G2 and due to (fT^ for positive Gi as well. 

Let us comment an essential point, connected with very small value of parameter uq. 
Note, that for m = solution satisfying all conditions excluding (fTTj) is given by the following 
expression 

20F °n ' ' 2' 2' ' ' 2/ ' 

with ratio 

r{u)u=o = G1/G2 = -3 +V12 = 0.464102. (14) 

At the same time value (fT^ . which follows from ()11|) is equal to 0.461538. For m — this 
ratio x{u) tends to value (fT^ from below. Because of these two numbers being quite close the 
intersection of curve x{u) with (fTT|) occurs at so small u. In case we use in expression flll|) the 
most leading terms in 1/A^, the ratio is 1/2 and there is no intersection of x{u) with this value 
at all, i.e. the solution does not exist. So a value uq is sharply dependent on approximations 
and quite small change of coefficients in expression pi|) influences it strongly. We take into 
account these considerations while commenting values of current quark mass in what follows. 

Solution Fi{p, —p,0,0) can be extended for more general kinematics. Namely, let us 
consider form-factor Fi{p, —p, q, —q). Assuming factorization property and bearing in mind 
the previous note on small contribution of Mq the form-factor reads as follows 



-p, g, -q) = ^Gtt{z\l, 1, 1, 0, 0, - ^) G«(/|l, ^, 1, 0, 0, - ^ 



^ = -^; ^=-^; x=p- y = q. (15) 

A study of form-factor F2{p'^), which enters into four-fermion vector and pseudo- vector terms 
in leads to result, that conditions of an existence of the corresponding solution does not 
provide additional restrictions for parameters being introduced above. Vector form-factors 
will be considered in details elsewhere. The explicit form of F2{p^) do not influence results 
of the present work. 



3 Wave function for scalar and pseudo-scalar excita- 
tions 

Now with the non-trivial solution of the compensation equation we arrive at an effective 
theory in which there are already no undesirable four-fermion terms in free Lagrangian (j2)) 
while they are evidently present in interaction Lagrangian Q. Indeed four-fermion terms 
in these two parts of the full Lagrangian differ in sign and the existence of the non-trivial solu- 
tion of compensation equation for Lagrangian (j21) means non-existence of the would be anal- 
ogous equation, formulated for signs of four-fermion terms in interaction Lagrangian 0. 
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In other words the fact, that sum of a series ^ G^a^ = 0, by no means leads to a conclusion, 
that sum of the same series with G ^ — G vanishes as well. 

So provided the non-trivial solution is realized the compensated terms go out from La- 
grangian Q and we obtain the following Lagrangian 

2 

Thus, bound state problems in the present approach are formulated starting from Lagrangian 

(uni). 

In case of realization of the nontrivial solution, let us consider zero-mass scalar and 
pseudo-scalar excitations . Let us write down Bethe-Salpeter equation for scalar zero-mass 
channel in the same approximation as was used in equation (jH)) 



ipT^-f^tp V'^SmV' + ^ AsTa^V'^^SsTm^) • (16) 



^, , GiN f°^^.^^ {Gl + QGiG2)N ( I r.^ o 2^ t / a ^ 

2 PX PX PX POO 

+ / (y) dy + '^og X / (y) dy + X log X / '^{y)dy+ / y log y (y) dy + 



+ X £ (log y + l)^{y) dy + f ^ ^ (^"^^ " ^ ^^^^ ~ 

-y y^{y)dy - log A^i^ J y^{y)dy + xj ^'(y)rfyjj; 

The corresponding differential equation for '^{x) is almost the same, as the previous one (|7j) 
with one essential difference. Namely the sign afore j3 is opposite. 



2 -1 + VI - 64u , -1 - VI - 64m /5m^ 
z = x- a = ; h= ; u=^. (18) 

Boundary conditions read 

PCO POO POO 

/ ^iy)dy = 0; / y^iy)dy = 0; / y^^iy)dy = 0. (19) 
J fi J fi J ^ 

Now we have four independent solutions decreasing at infinity, which we use for general 
solution 

^(^) = C*Gf,[^\l^ i, 0, i, a, h) + CrG^°(^|l, i, i, 0, a, fo) + 

+ ClGf,[^\l^ a, 6, i, i, O) + ClGf,[^\\^ a, 6, 1, i, o) . (20) 



Boundary conditions (fT^ allows to express C*, i = 1,2,3 in terms of C* and condition of 
wave function to be unity on the mass shell ^>{m?) = 1 fixes C* as well. The normalization 
condition of the Bethe-Salpeter wave function fixes constant of interaction of the bound state 
with quark-anti-quark pair 

5( (</)^^ + z7ra^75raV) ■ (21) 



Namely this normalization condition, which corresponds to a correct coefficient afore a kinetic 
part of scalar and pseudo-scalar particles, gives 

The loop integral of the fourth order in interaction (j?H) gives four-fold interaction 

a(^1^; (23) 

where 

rm±. (24, 

The result, that equation (fTTj) has unique solution, which satisfies all boundary condi- 
tions, is the confirmation of the evident fact, that in the same approximation, in which the 
compensation equation has a non-trivial solution, fields (f) and tTq are to have zero masses 
in accordance with Bogolubov-Goldstone theorem [21 El El- However an account of chromo- 
dynamic interaction leads to additional contribution to the masses. Let us calculate a mass 
correction term due to QCD interaction. For the purpose let us take into account terms of 
the first order in P^, where P is the momentum of a scalar (and pseudo-scalar) meson and 
one-loop QCD term. We have 

^^^^ - ^ J + 32^^^ 

/ n2i {P-Q? 3, ,^\/ 3P2 (qPf^^^/(q^)dq 
2A^+(p-,)^log^-5(p-#) (1-^ + 5^)^ + 

+ ^ fim^^. (25) 

4vr4 y g (g -p) 

In the course of QCD term calculation we use transverse Landau gauge ^ . Let us multiply 
equation by '^{jP')/p'^ at P = and integrate by p. Due to equation (fTTj) be satisfied 
we have 

P2 [^>{q^fdq gl f m{p^)dp f m{q^)dq 

+ 731 / 3 / ..J. = ; (26) 



2 J (g2)2 ' 4 7r4 7 J q^{q-pY 



^In the approximation used the transverse gauge leads to absence of renormalization of both vertex and 
spinor field 
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After angle integration we get 



A ^ 7„ ' 64 ' 64 ■ ^ ^ 

The integral entering into ^TJ\ looks like 

+ Q'^06(^lf,a^6^0,|,l) + Q'^06(^ll,a',fe',0,|,f) ) ! (28) 

/ u 

1 - VI - 64m 1 + VI - 64m 

4 ' 4 

Relation (P7j) gives us mass of scalar (and pseudo-scalar) in terms of average low-energy 
QCD constant a^. There are a number of considerations concerning possible values of the 
parameter. For example, in approach [Hj with low-energy freezing of the strong constant 
value (is = 0.414 is quoted. Smooth matching of perturbative and non-perturbative regions 
in QCD gives value dg = 0.4354 fTUJ. There is also a sum rule definition of this parameter 
from experimental data dg = 0.47 ± 0.07 there are also considerations on behalf of 

larger values of ds 113j . Taking into account these remarks we consider range of values 
of the effective low-energy constant ds = 0.4 — 0.75. In other words we can define the same 
quantity as value of running strong constant at a characteristic momentum, e.g. 600 MeV. 
Different variants of low energy behaviour of as lead to the same interval. Thus we have 

(29, 

where Ji is defined in fl22|l . Due to both integrals being positive we evidently have tachyon 
mass and so a scalar condensate appears in the minimum of effective potential 

* 2 24 

For the present approach it is highly important, that the account of the QCD interaction 
leads to tachyon mass of a scalar and pseudo-scalars. The appearance of tachyons and thus 
the appearance of the scalar condensate in the minimum of the effective potential results in 
stability of the non-trivial solution and consequently one may conclude, that Lagrangian p6|) 
with the conditions being obtained are valid. 

Bearing in mind definitions of quantities entering in the effective potential, we obtain for 
value of the scalar condensate rj 

— 6 Tfl^ 

r/2 = ■ T] =<(f)> . (30) 
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Mass of after symmetry breaking is iy2(^^^mfy, and vr mass equals to zero in the present 
approximation. The constituent quark mass ruq is defined by relation 



mg-mo=gv = 9\l ^^V^' 

Integrals in this as well as in other relations of this section including boundary condi- 
tions are calculated with the lower limits defined not by the current mass mo, but by the 
constituent mass m. Now parameter u is defined by relation 

Thus after the appearance of the scalar condensate the quark propagator takes form 

G{v) = ^^^^ [ ^^^^ ; = (m - mo) ^(p') + mo . (33) 

Relations (jSH EH EHI) give for a, = 0.46 

M = 0.001; C* = -0.919; C* = -0.0255; C* = - 1.895 ; 

C; = 0.228 ; Ji = 1.285 ; h = 0.801 ; h = 0.302 . (34) 

We take cxg and current mass mo as the initial parameters and express in their terms 
all other parameters including vr-decay constant /,r- We use for the latter the Goldberger 
- Treiman relation. In the framework of the present approach we obtain the relation by 
considering the expression for transition 7r+ fi'^ u^. We have 

A = f4 / ({m-mo)^{yf + mo^iy))— = 

'^'^ J vn? ' y 

= |^((m-mo)/i+ mo/rj ; h = j • (35) 

Provided either mo = or Ji = J7 we get with account of normalization condition ()22p just 
the original Goldberger - Treiman relation m = g f.^. We use full relation . However let 
us note, that values of the two integrals are close Ji ~ J7 and the simple original relation 
works with sufficient accuracy. 

At this place it is worth-while to explain that we have a solution for any value of mo and 
the problem is, which value is to be chosen. While mo is not measured directly, the accuracy 
of its definition is not very high. For comparison with experimental data we propose to fix 
value of /tt = 93 MeV, which is quite certain, while other parameters including mo being 
subjects for calculation. 

As a result we obtain e.g. for = 0.46 {u = 0.001) 

9 = ^-'^-^Gi= (240 5 MeV)2 ' = 298.5 MeV; mo = 19.8 MeV. (36) 
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Mass of a is defined by above results IH^ 



= 2(m-mo) W-^ = 440.1MeV. (37) 
V h 

At this stage mass of vr-meson is zero due to appearance of vacuum average of scalar field 0. 
Coupling constant ga-Ki, in interaction 

(7,^^(T7r'^7r^ (38) 

is calculated with the use of wave function (pn|) with parameters (1341 IH^ . The triangle 
diagram gives 

= ^ / ^^^((m-mo)vl/(i/)+mo) = 8.04m. (39) 

According to ipUj) the a width reads 

3 0^ 

= ^'^^^ v^m2 - 4m2 = 605.8 MeV; = 139 MeV. (40) 

16 vrm^ 

Values (j?71 HOI) do not contradict the existing data [H], which, as a matter of fact, are 
characterized by a wide spread of results. 

The results on the wave function of a scalar particle and on the form-factor of the effective 
four-fermion interaction allow to estimate the quark condensate in the next approximation. 
Indeed vertex aip ip form- factor is just g'^{q'^) . The non-perturbative part of the quark prop- 
agator reads gri'^{q^) = (m — mo)\E'(g^) . According to definition of the quark condensate 
we have 

<gg> = rW)dq\ (41) 

Here the integral equals zero due to boundary conditions (jl9p . However one may try to 
calculate subsequent approximations. Let us calculate loop corrections: firstly with vertex 
corresponding to four-fermion interaction in Lagrangian ()16|) and secondly with gluon ex- 
change in transverse gauge (see Fig. 3). In loop diagrams with the form-factor Fi{p) we sum 
up infinite loop chain (the senior orders in 1/A^), that gives factor (l — B) in the denominator 
of the first term in the following expression 

12y/3{m - mo) 3as{m - mo) f ^ 16 vr^ 

„ 16 riFiizyfdz 



VU Ju 
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where function "^liz) is defined by relation ()28|) . We have e.g. again for u = 0.001, = 0.46 



that is sufficiently lower than the traditional value (240 MeV)'^. Nevertheless we draw at- 
tention to the correct sign of the condensate. It is connected with the sign afore Gi in the 
interaction Lagrangian ()16p .which is defined by conditions of the existence of the non-trivial 
solution of the compensation equation. 

4 Discussion of results 

The low-energy hadron parameters being calculated for three values of are presented 
in Table 1. The calculations needs numerical integrations of expressions containing Meijer 
functions with infinite upper limits. Therefore the presented data may contain uncertainties 
(not more than 3%). Note that integrals /4, Jg, It are evaluated analytically, and all the 
rest ones are calculated numerically. Let us recall, that only two parameters are in our 
disposal - the average strong coupling constant and dimensional constant of vr-decay /^r. 
We use of course the Goldberger-Treiman relation (|35p. which is obtained in the framework 
of the present approach. No other outside information is used. In particular while calculat- 
ing (j-meson width we substitute pion mass which was obtained here in the corresponding 
approximation and was quoted in the Table. We see from Table 1, that the approach gives 
reasonable correspondence to the existing facts and data. We present information from ex- 
periments and from phenomenological considerations in the last right-hand column of the 
table. 

Value of current mass mo is essentially larger than customary values. This parameter 
is defined by value uq fll2|) . which is sharply dependent on corrections of the next orders of 
1/A^ expansion and on other details of the next approximations. For example, by a slight 
change in coefficient afore the second term in brackets in expression (II Ij) one may obtain 
for mo value being two-three times smaller than that presented in Table 1. At this all other 
parameters practically do not change. Taking into account this fact we do not consider the 
deviation of mo as critical for the approach. In comparison to the traditional value the 
modulus of quark condensate is essentially smaller. While its calculation we have seen that 
in the leading approximation this parameter is zero and the values quoted correspond to loop 
corrections. It is possible, that there are other contributions to the quark condensate, which 
give essential changes. For example contribution of s-quark loop to the first term in ()42|) 
evidently enlarge it. This problem deserves a special study. However just the presented values 
give rise to an interesting observation. It comes out, that Gell-Mann - Renner relation O^j 
agrees sufficiently well with the parameters of Table 1. Indeed In the lowest line of the table 
for pion mass we write down its values calculated by the relation 



< qq > 



(167 Me V) 



3 . 



(43) 



2 



2 mo 



m: 



7r 



J TT 



< qq > . 
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Values of for mo and < qq > from the table are quite satisfactory in the range of a 
reasonable accuracy. Note also that values of the quark condensate being lower than the 
traditional one are discussed in the literature (see e.g. jtt^). Lattice measurements [T7j 
give rather wide interval of renormalization invariant value of the condensate — < qq > 
= ((206 ± 44 ± 8 ± 5) MeV)^, in which also the values presented in Table 1 also enter. 
Parameters of the a-meson especially its width noticeably depend on a choice of value as- 
Unfortunately the spread of data [H] covers all presented results. Provided one choose results 
of one work then it is possible to do some conclusion on preferable value of a^- For example 
in paper [TH] the processing of the set of data leads to result = (470 ± 30) MeV, = 
(590 ± 40) MeV. Such values agree well with our results for ~ 0.5. At that the pion mass 
and the constituent quark mass are quite satisfactory. 

So we may state a quantitative agreement for constituent quark mass and for parameters 
of mesons vr and a, whereas for current mass and for the quark condensate we may declare at 
least a qualitative agreement. One hardly could expect more from the first non-perturbative 
approximation with only two parameters, one of which /^^ is strictly fixed. Therefore one 
have to admit the description of the low-energy data by the approach being satisfactory. 

To conclude let us emphasize that the aim of the work is achieved. We have begun with 
the demonstration of the non-trivial solution of the compensation equation. The existence 
of scalar and pseudo-scalar excitations (mesons) in the same approximation is a consequence 
of its existence. The account of QCD interaction leads to the shift of their masses squared 
to the negative region, i.e. to the appearance of tachyons, which are necessary for scalar 
condensate to arise. As a result we obtain the standard scheme of chiral symmetry breaking 
with massive scalar and massless pion. Subsequent approximations of the approach are 
related to values of the quark condensate and the pion mass. 

We have shown that the application of the method of work T], which is based on Bo- 
golubov quasi-averages approach, to the low-energy region of hadron physics leads to quite 
reasonable results. From this we would make two essential conclusions. 

Firstly, a subsequent development of the present approach to the hadron physics quite 
deserves attention. In particular it is advisable to apply the approach to calculation of 
parameters of vector mesons p, u, Ai, to consider hadrons containing s-quark, to take into 
account the vr — Ai-mixing, to introduce diquarks etc.. In view of methods it would be 
desirable to improve a procedure of taking into account of particle masses. These problems 
comprise subjects for a forthcoming studies. 

Secondly, the positive result of applicability test with Nambu - Jona-Lazinio model, 
being taken as an example, allows to hope for successful application of the approach to other 
problems. In particular we mean the problem of a dynamical breaking of the electroweak 
symmetry. A qualitative discussion of possible variants in this region is presented e.g. in 
works [in], 120] • 

The author is grateful to M.K. Volkov and S.B. Gerasimov for valuable discussions and 
to I.V. Zaitsev for some numerical calculations. 

The work is partially supported by grant "Universities of Russia" UR.02.02.503. 
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A Appendix 



We present below a set of formulas for angular integrals in four-dimensional Euclidean space, 
which are used in the work. 



X = , y = q^] 



TT^ / ydyF{y)(d{x-y)- + d{y-x)-)- 



d^qF{q^) 

{p-qf " Jo ^-^^^ ^^^V"'"' ^'x ' "^'y. 

'himM = ryayF(y)(e(.-y)'l^e(y-.)^ 

[p-q) 2 Jo \ X yJ 

dSF{q^){pqf roo , ,„.2 



{p-qf 
d^qF[(^)[j)qf it 
{p - qf 

d'qF{q'){p,p-q) 2 T ^ W( W 
dWq^)iq,q-p)_^, r...W(.-.)-^ 



— ydyF{y) [e{x ' y)\^ + v) + ^iv - ^) + — 
^ ydyF{y) {9{x -y){^ + ^y') + Oiy - x) (2x' + ^ 



d^qFmp,p-q)ipq) vr^ T ^,^^(^) ^,(, _ ^) % _ _ 



((p-g)2)2 4 7o ^ X y 

j d^qFiq') log (g - p)^ = vr' ^ (^(x - Z/) + log x) + 

+ ^(y-x)(^iog y + ; 

2 

+ ^^(y-x)(-3x + ^)); (44) 
d^qF{q^){pqf log (q - pf = — J ydyF{y) [e{x ^ + ^2/ log + 



+9{y - x) \xy\og 2/ + ^ 
j d'qF{q^){pqf log (g - P)' = y ydyF{y) {e{x - + ^ - 2y2a;) + 

+ eiy - x)i^-2x^y + — + —jj . 
Integrals containing log {q — pY are evaluated using formulas from book [21J. 
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Tabic 1. 

Parameters of the low-energy hadron physics for three values 
of effective coupling constant ocs- 



Ots 


0.46 


0.61 


0.73 


exp/phen 


AMeV 


93 


93 


93 


input 


9 


3.20 


3.44 


3.59 




mo MeV 


19.8 


17.9 


16.9 


5-10 


mMeV 


298.5 


321.3 


335.5 


300 - 350 


Q-^''^ MeV 


240.5 


217.6 


205.4 




- <qq >^/^MeV 















167.0 


155.4 


148.7 


220 - 240 


m^r MeV 















146.0 


124.6 


113.3 


136 - 140 




8.04 


8.27 


8.40 




rria MeV 


440.1 


469.1 


486.3 


400 - 1200 


r„ MeV 


581.1 


761.1 


862.5 


600 1000 


u 


0.001 


0.002 


0.003 
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Figure captions 



Fig. 1. Diagram representation of tlie compensation equation. Black spot corresponds to 
four-fermion vertex with a form-factor. 

Fig. 2. Diagram representation of Bethe-Salpeter equation for scalar bound state, the 
later corresponding to a double line. 

Fig. 3. Loop correction to the non-perturbative part of the quark condensate. The full hne 
corresponds to the quark propagator with non-perturbative mass operator (m — mo)^(p^). 
The dotted line represents a gluon. 
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Fig. 



